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Abstract 

Equation of long-range particle drift and diffusion on three-dimensional physical lat¬ 
tice is suggested. This equation can be considered as a lattice analog of space-fractional 
Fokker-Planck equation for continuum. The lattice approach gives a possible microstruc- 
tural basis for anomalous diffusion in media that are characterized by the non-locality of 
power-law type. In continuum limit the suggested three-dimensional lattice Fokker-Planck 
equations give fractional Fokker-Planck equations for continuous media with power-law 
non-locality that is described by derivatives of non-integer orders. The consistent deriva¬ 
tion of the fractional Fokker-Planck equation is proposed as a new basis to describe space- 
fractional diffusion processes. 

PACS: 05.20.-y; 52.65.Ff; 45.10.Hj; 61.50.Ah 


1 Introduction 


Fokker-Planck equations are usually used to describe the Brownian motion of particles [T]. 
These equations describe the change of probability of a random function in space and time in 
diffusion processes. The Fokker-Planck equation is usually the second-order partial differen¬ 
tial equation of parabolic type. In many studies of diffusion processes in complex media, the 
usual second-order Fokker-Planck equation may not be adequate. In particular, the probability 
density may have a thicker tail than the Gaussian probability density, and correspondent corre¬ 
lation functions may decay to zero much slower than the functions for usual diffusion processes 
resulting in long-range dependence. This phenomenon is known as anomalous diffusion [2l[3lll]. 
Anomalous diffusion processes can be characterized by a power-law mean squared displacement 
of the form [21 [3l 0] 


{x^{t)) 


2K{a)r 
r(a + i)’ 


( 1 ) 


where r( 2 ;) is the Gamma function, a is the anomalous diffusion exponent, and K{a) is the 
anomalous diffusion constant. In equation ([1]), we use the second moment that is dehned in 
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terms of the ensemble average. Depending on the value of a, we usually distinguish sub¬ 
diffusion for 0 < a < 1 or super-diffusion for a > 1. There are two limit cases such as 
the normal diffusion (a = 1) and the ballistic motion {a = 2). One of possible approach to 
describe the anomalous diffusion is based on the continuous time random walk models [5] in 
which the particles are considered as random walkers with step lengths r and waiting times 
t. An important role play the anomalous diffusion processes with the Poissonian waiting time 
and the Levy distribution for the jump length. The Levy flights [2] are random walks in which 
the step-lengths (long jumps) have a probability distribution that is heavy-tailed. The Levy 
motion can be described by generalized diffusion equation with space derivatives of non-integer 
orders /r, |1]. The fractional moment of order S for Levy flights has the form (|a:(f)|^) 
where 0 < ^ < p < 2. 

Derivatives of non-integer orders [SI El El El CHI IHl [121 US] play an important role in describing 
particle transport in anomalous diffusion jH |ITl |15l ESI EZl Ell EH] and have a wide application in 
various areas of physics (see for example [201 IHIl IHHl |23l [Ml IHSl IHSl 123128]). Various approaches 
lead to different types of space-time fractional Fokker-Planck equations. Usually the space- 
fractional Fokker-Planck equations are obtained from the second-order differential equations 
by replacing the hrst-order and second-order space derivatives by fractional-order derivatives 
[7]. Fractional Fokker-Planck equations with coordinate derivatives of non-integer order have 
been suggested in [29]. The solutions and properties of these equations are described in [30l 
ITS] . The Fokker-Planck equation with fractional coordinate derivatives was also considered in 
[SB El E21 EH ES ES]. It should be noted that the fractional Fokker-Planck equations can 
be derived from the probabilistic continuous time random walk [SEl EB ES]- In this paper 
we propose a consistent derivation of the space-fractional Fokker-Planck equation based on 
lattice model with long-range drift and diffusion that is considered as a microstructural basis 
to describe fractional diffusion processes in continua. 

A discrete lattice version of the Fokker-Planck equation in analogy with the lattice-Boltzmann 
models has been suggested in [IHl SB S2] • These models are used to solve the equations of hy¬ 
drodynamics and cavity flow simulations [13] . The lattice Fokker-Planck equation is applied to 
the study of electro-rheological transport of one-dimensional charged fluid [H] and it is used 
in phase-space description of inertial polymer dynamics [15]. All these lattice Fokker-Planck 
equations are based on the lattice Boltzmann discretization approach. 

In this paper, we propose a lattice equation for probability density of particle in unbounded 
homogeneous three-dimensional lattice with long-range drift and diffusion to n-site from all 
other m-sites (m ^ n). We prove that continuous limit for the suggested lattice Fokker-Planck 
equation gives the space-fractional Fokker-Planck equation for non-local continuum. The frac¬ 
tional differential equation for continuum contains generalized conjugate Riesz derivatives on 
non-integer orders. 

Continuum mechanics [IB] can be considered as a continuous limit of lattice dynamics 
[IBllHlllgllso], where the length-scales of a continuum element are much larger than the 
distances between the lattice particles. The hrst self-consistent derivation of the Fokker-Planck 
equation based on the microscopic dynamics for classical and quantum systems was obtained by 
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Bogolyubov and Krylov (SUES]. Long-range interactions are important for different problems 
in statistical mechanics [531 El EH], kinetic theory and noneqnilibrinm statistical mechanics 
[561EZ], theory of non-eqnilibrium phase transitions jSHlEH]. As it was shown in [6lEn] (see also 
[HEHM] and |65l El El ESI Ell ED] ), the continunm equations with fractional derivatives can 
be directly connected to lattice models with long-range properties. A connection between the 
dynamics of lattice system of particles with long-range properties and the fractional continuum 
equations are proved by using the transform operation [HI ED]- The papers [HI ED] deal 
with the one-dimensional lattice models and the correspondent one-dimensional continuum 
equations. In this paper, we suggest three-dimensional lattice models for space fractional 
diffusion processes. We propose a general form of three-dimensional lattice Fokker-Planck 
equation, which leads to continuum fractional Fokker-Planck equation with space derivatives of 
non-integer orders by continuous limit. The suggested approach to derive the space fractional 
Fokker-Planck equations can serve as a microstructural basis to describe the spatial-fractional 
diffusion processes. 

2 Lattice with long-range drift and diffusion 

The lattice is characterized by space periodicity. In an unbounded lattice we can dehne three 
non-coplanar vectors ai, a 2 , as, that are the shortest vectors by which a lattice can be dis¬ 
placed and be brought back into itself. All space lattice sites can be dehned by the vector 
n = (tt-i, n2, n-s), where n* are integer. For simplihcation, we consider a lattice with mutually 
perpendicular primitive lattice vectors ai, a 2 , a 3 . We choose directions of the axes of the Carte¬ 
sian coordinate system coincide with the vector a*. Then a* = a* e*, where Oj = |aj| and e* are 
the basis vectors of the Cartesian coordinate system. This simplihcation means that the lattice 
is a primitive orthorhombic Bravais lattice with long-range drift and diffusion of particles. 

If we choose the coordinate origin at one of the sites, then the position vector of an arbitrary 
lattice site with n = (ni,n 2 ,n 3 ) is written r(n) = niai -|- n 2 a 2 -|- n 3 a 3 . The lattice sites are 
numbered by n, so that the vector n can be considered as a number vector of the corresponding 
particle. We assume that the positions of particles coincide with the lattice sites r(n). The 
probability density for lattice site will be denoted by /(n, t) = f{ni,n 2 ,n 3 ,t), where the site is 
dehned by the vector n = (rii, 77 - 2 , n 3 ). The function /(n, t) satishes the conditions 

-j-OO -hcso -1-00 

E E E /(m, 772 , 773 , t) = 1, /(77i, 772 , ^ 3 , t) > 0 (2) 

72,1= —00 712 = —00 n3 = —00 


for all f G M. 

The equation for probability density of particle in unbounded homogeneous lattice is 


dt 


3 3 

= ^ ^ ^ (n - m)/(m, f), (3) 

i=l rrii^ni i^j=lmi^ni rrij^nj 
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where /(n, t) is the probability density function to find the test particle at site n at time t. 
The italics i,j G {1;2;3} are the coordinate indices, gi and gij are lattice coupling constants. 
The coefficients — m) and ^.(n — m) describe the particle drift and diffusion on 

the lattice, and it can be called the drift and diffusion kernels for lattice step length n — m. 
These kernels describe the long-range drift and diffusion to n-site from all other m-sites. The 
parameters ctj and f3j in the kernels are positive real numbers that characterize how quickly the 
intensity of the drift and diffusion processes in the lattice decrease with increasing the value 
n — m. These parameters also can be considered as degrees of the power law of lattice spatial 
dispersion [65l [68] that is described by non-integer power of the wave vector components. 

Equation Q describes fractional diffusion processes on the physical lattices, where long- 
range jumps can be realized. The Levy motion (flights) for these lattices can be described by 
the lattice Fokker-Planck equation (|3|), which is considered as a lattice analog of the fractional 
diffusion processes with the Poissonian waiting time and the Levy distribution for the jump 
length [Ij. 

For simplihcation, we consider the kernels in the form 

(n - m) = {Ui -rrii), (n - m) = - rrii) Kp.{rij -mj), (4) 

where i,j = 1,2,3. The kernels ~ "^t), where i = 1,2,3, describe long-range jumps in 

the direction sg with lattice step length — rrii in the lattice. The correspondent terms with 
kernels K^X'^i ~ ^i) can be considered as lattice analogs of fractional derivatives of order aj 
with respect to coordinate Xi = (r, a*). We will consider even and odd types of the kernels 
KaX'^^i ~ nr*), i = 1,2, 3, that will be denoted by — nrj) and K~Xni — rrii) respectively. 

We assume that the kernels KX{n) satisfy the following conditions: 

1) The kernels iP^(n) are real-valued functions of integer variable n G Z. The kernels KX{n) 
and K~{n) are even and odd functions such that 

K(-n) =+Kt(n), K-{-n) = -K-{n) ( 5 ) 

hold for all n E 

2) The kernels belong to the Hilbert space of square-summable sequences, 

OO 

'^\KX{n)\‘^ < oo. (6) 

n=l 

3) The Fourier series transforms K^ik) of the kernels in the form 

+ 050 050 

X] e““”‘At(") = 25;A^(n)cos(fa), (7) 

n=—oo n=l 

+ CO 050 

K-{k) = ^ e-““A--(n) = -2 J ^ K-{n) sm(te) (8) 

n=—oo n=l 
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satisfy the conditions 


K*{k) = \k\-+o{\kn, (k^o), 


(9) 


and 

-^« (^) = * sgn(A;) |/c|“ + o(|/c|"), (A;-)■ 0) (10) 

respectively. Here the little-o notation o(|A;|") means the terms that inclnde higher powers of 
\k\ than |/c|". The snggested forms ([9]) and flTOl) of the Fonrier series transforms of the kernels 
K^{n) mean that we consider lattices with weak spatial dispersion [65]. The conditions (I9|) 
and ffTOj) allow us to consider a wide class of kernels to describe the long-range lattice drift and 
diffusion. 

In general, type of dependence of the function on the wave-vector k is defined by 

the type of the spatial dispersion in the lattice [651 EH]- For a wide class of processes in the 
lattice, the wavelength A holds the relation Oq/A ~ ka^ -C 1, where Oq is the characteristic size 
of the lattice distance such that oq = max{|ai|, |a 2 |, |a 3 |}. In the case kao <C 1, where oq, the 
spatial dispersion of the lattice is weak. To describe lattices with such property it is enough 
to know the dependence of the function only for small values A;, and we can replace 

this function by the Taylor’s polynomial series. The weak spatial dispersion of the lattices with 
power-law type of spatial dispersion cannot be described by the usual Taylor approximation. In 
this case, we should use a fractional Taylor series [65116H] . The fractional Taylor series is more 
adequate for approximation of non-integer power-law functions. For example, the usual Taylor 
series for the power-law function K^{k) = A:" has the infinite by many terms for non-integer 

a. The fractional Taylor series of order a has a hnite number of terms for this function, and 
the fractional Taylor’s approximation is exact. We can use the fractional Taylor’s series in the 
Riemann-Liouville form (see Chapter 1. Section 2.6 m) that can be represented as 

KUk,) = b,{a)\kr + o{\k,n, ( 11 ) 


where 


bj (a) 


F(cr -|- 1) 


( 12 ) 


and qD'^ is the Riemann-Liouville fractional derivative [7] of order 0 < a < 1 with respect to 
k. This derivative is defined by 




m 



(k), 


( 13 ) 


where is the left-sided Riemann-Liouville fractional integral of order a > 0 with respect to 
A) of the form 

= («=> 0 ). ( 14 ) 

Using the approximation flTT|) . we neglect a frequency dispersion for simplification, i.e., the 
parameters bj{a) do not depend on the frequency u. In suggested lattice models, we define the 
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kernels such that the constants gi and gij include the factor bj{a), and as a result the conditions 
(E]) and flTOl) hold. 


For simplihcation, we can consider the lattice kernels that are dehned by the explicit ex¬ 
pressions in the form 

K+ik) = l^r, k-ik) =tsgn{k) \kr. (15) 

In this case, the inverse relations to the dehnitions of k^{k) by equations ([7]) and dS]) have the 


forms 


ik — ~ cos(n k) dk, {k — -/ sin(n k) dk. 


TT 


TT 


(16) 


For non-integer real values of the parameter a, the expressions for the kernels K^{n — m) are 


K^{n — m) = 


TT 


Oi \ 


1-^2 


a + 1 1 a-|-3 7r^(n — m)^ 


7 r"+^ (n-m) ^ 

K (n - m) =- 


2 2 2 4 

a+ 2 3 a-1-4 7r^(n — m)^ 


a > — 1 , 


a > - 2 , 


( 17 ) 


( 18 ) 


where 1 F 2 is the Gauss hypergeometric function (see Chapter II in [72]). Note that expressions 
can be used not only for a > 0, but also for some negative values of a. 

To visualize the properties of the kernels ffTTll and ffT^ . we give the plots of the functions 


F+{,x,y) 


TT^ fy + lly + ?> 

5TT (“’ 2- ~~j' 


( 19 ) 


where 


F_{x,y) 


X fy + ‘2 3 y + ^ x‘^\ 

y + 2 ' ' V^’ 2 ’ ^) ’ 


K^{n — m) = F±{n — m, a). 


( 20 ) 

( 21 ) 


We present the plots of the function flT^ by Figures 1, 3, 5, and the plots of fl20|) by Figures 2, 
4, 6 for the same ranges of x and y > 0. 

Let us note some qualitative properties that can be seen from Figures 1-6. We should note 
that the functions (IT^ and (j^ represent the kernels with (1T7|1 and (1T8|1 that describe the 
long-range drift and diffusion to n-site from all other m-sites, where m G N. Oscillations tell 
us that the inflow and outflow of probability periodically change each other, when the distance 
X = n — m between sites increases. The negative values of F±{x, a) can be interpreted as the 
probability flux from the site, and the positive values of F±{x,a) can be interpreted as the 
flux to the site. Maximums and minimums of F^{x, a) characterize an amplitude of oscillation 
of the probability flux from the site and into the site. The amplitudes as functions of the 
parameter a are increasing functions for a hxed value x = n. Plot of the functions fll9p and 
fl20|) with a = 1.5 and a = 1 for the range a; G [0,7] are presented by Figures 7 and 8, where the 
graphics of functions with a = 1.5 have larger amplitudes than the graphics of the functions 
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with a = 1. The amplitudes as functions of the values n are decreasing functions for a hxed 
value a, and this decreasing has a power-law form. During the transition from non-local to 
local case for the functions F±{x, a), a sharp jump does not occur. We can only state that the 
fractional power-law decreasing is transformed into the decreasing of integer power form. 

It should be noted that the kernels (n) give the local operators for continuum limit for 
even a only, and K~{n) give the local operators for odd a only. The kernels K^{n) for integer 
values of a (see also Sec.2.5.3.5 in m) can be represented by the equations 


Kin) 


[(a p/2] ^ (-l)[(“+b/2]«!(2[(a-t-l)/2] 

(a — 2n — l)\ fi 2 k+ 2 7 r? 7 ,"+^ 

k=0 ^ ' 


a) 


( 22 ) 


and 


[a/2] 




{_l)n+k+i a\ TT^^-^k-i i (-1)a! (2[a/2] - « + 1) 


k=0 


{a — 2n)! 


n 


2k+2 


Tin 


Qf+l 


(23) 


where [z] is the integer part of the value z, and 2[(a -|- l)/2] — a = 1 for odd n, and 2[(a 
l)/2] — a = 0 for even n. We can give examples of kernel with some integer a. Using equation 
fl2^ or direct integration ffT6|) for a G {1; 2; 3}, we give K'^{n) in the form 


K*(n) = 


1 - (- 1 )' 




A'+(n) = 


2 (-l)' 




Kt(n) = 


3 ir(-l)" 6(l-(-l)") 




71 


(24) 


h7(n) = 


- 1 )' 


n 


A'-(„) = (zi):i+!(iUun, 


n 


TTn-^ 


K7(n) = 


(-l)”ir2 6(-l)' 


n 




, (25) 


where (1 — (—1)"') = 2 for odd n, and ((—I)” — 1) = 0 for even n. 

For simplihcation of the form of lattice equation, we use the lattice operators ["] such 
that the action of these operators on the lattice probability density /(m, f) is 


a 
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+00 

Kini-rrii) 


mi=—oo 

rrii^ni 


(i = 1,2,3). 


(26) 


The values i = 1,2,3 specify one of the three variables ui, U 2 , of the lattice site that are 
similar to Xi of the space If Oj = 1, then K+ is a nonlocal operator, and if cxj = 2, then 

Mr are nonlocal operators also. Note that the operators K+ ["*] for odd integer values of ai 
and K“ [“*] for even integer values of ai are nonlocal. For example, the operators [^] and 
M~ [^] cannot be considered as a local operators of integer orders. 

We also can consider combinations of the lattice operators 



aiPj 

= 

'ay 


■ 5^ 

_ 1 


. * i . 


- i . 


IJ J 


( 27 ) 
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where i, j take values from the set {1;2;3}. The action of the operator fITrj) on the lattice 
probability density /(m, f) is 




ai(3j 

i j 




+ 00 

E 


+00 


mi=—oo mj=—oo 
rrii^rii rrij^rij 




(28) 


This is the mixed lattice operators. 

Using the lattice operators fl26l) and fl28|l . the equation for probability density ([3]) takes the 
form 


dt 


i=\ 


*P=i 


±,± 


oti (5 j 

i j 




(29) 


This is the three-dimensional lattice Fokker-Planck equation in the operator form to describe 
fractional diffusion and drift with the lattice jump length (n — m). 

To describe the long-range drift and diffusion for the lattice with memory, we can use the 
equation 


a/(n,f) 

dt 


^9i R 

i=l 


zb 


(Xi 


L I 




9ij IK 


±,± 


13j 

* j 


RL 

0 


D] ^/(m,f), 


(30) 


where is the Riemann-Liouville fractional derivative of order (1 — 7 ) with respect to 

time [7]. Note that the time-fractional derivatives is present only in the diffusion term. 

This fractional derivative describes the long-term memory of power-law type. Equation fl3UI) 
describes anomalous diffusion processes with the waiting time t and the lattice jump length 

(n — m). 

We can consider the time-fractional derivatives in the hrst and second terms of the 

right side of the lattice Fokker-Planck equation fl29|l . In this case the time-fractional lattice 
Fokker-Planck equation has the form 


.a, p 

— n -Lyx i-' T T? 


dt 


LFP 




where C 


q.,13 

LFP 


is the lattice Fokker-Planck operator 


_ 

^LFP ~ 


3 


'^9i 


■ay 

- i - 


3 


+ 9ij 

*j=i 


ai^j 

i j 


(31) 


(32) 


Equation (|3Til describes long-range diffusion and drift with power-law memory on orthorhombic 
Bravais lattices. 
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3 Continuum limit for lattice equations 

3.1 Continuum limit for lattice probability density 

In order to transform a lattice probability density /(n, t) into a probability density f{r,t) of 
continuum, we use approach suggested in [601161]. We propose to consider /(n, f) as Fourier 
series coefficients of some function /(k, f) for kj G [— /cjo/2, fcjo/2], where fcjo = 27r/aj. Then 
we use the continuous limit ko —)■ cxo to obtain /(k, f), and finally we apply the inverse Fourier 
integral transformation to obtain the probability density f{r,t). For clarity, we have presented 
the set of transformations of the probability density by Figure 9. 



Fourier series 


Inverse Fourier 

transform J^a 


integral transform 



Figure 9: Transformation of lattice probability density to continuum probability density. 


The transformation a lattice probability density into a continuum probability density is 
realized by a sequence of the following three steps: 

The first step is the Fourier series transform : /(n, t) —)■ J^A{/(n, t)} = /(k, t) that is 

defined by 

+ 0O 

/(k,()= 5^ (33) 

ni,n2,n3=—oo 
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where the inverse transformation is 


/ 1 \ /■+*^io/2 r+f^3o/'^ 

= IIt- / / <'*'3 /(k,() =Pi'{/(k.()}. (34) 

\ .-1 ^.iO I J—k-\r)/2 J —k9,c\/2 


f+fc3o/2 




—fclo/2 


—^ 30/2 


and r(n) = ^jo = We assnme that all lattice particles have the same 

inter-particle distance in the direction aj for simplification. 

The second step is the passage to the limit aj —)■ 0 {kjo —)■ oo) denoted by him : /(k, t) —?• 

Lim{/(k,f)} = f{k,t). The fnnction /(k,f) can be derived from f(k,t) in the limit a* —)■ 0. 
Note that /(k, f) is a Fonrier integral transform of the probability density f{r,t), and f{k,t) 
is a Fonrier series transform of f{n,t), where we nse 



considering r(n) = = 27r nj/fejoe^ —)■ r. 

The third step is the inverse Fourier integral transform : 

/(r, t) is defined by 

1 r+°° 3 

/("•*) = 7W / 

{27ry 

that corresponds to the transformation 

/ + 00 

-OO 

Note that the Fourier series transform equations fl3^ and fl3T|) in the limit aj — )■ 0 {kjo — )■ cxd ) 
give the Fourier integral transform equations fl36|) and fl35|) . where the sum is replaced by the 
integral. 

The lattice probability density /(n, f) is transformed by the combination o Limo 
into a probability density f{r,t) of continuum, 

o him o (^/(n, t) j = /(r, t). (37) 

The combination of the operations Lim, and allows us to map the lattice functions 
and operators into functions and operators for continuum. 


/(k,f) = 

ViKt)} (35) 

r,f)}. (36) 


3.2 Continuum limit of lattice operators 

Let us consider transformations of lattice operators (126]) and (127p into continuum operators. 
The transformations o Limit o map the lattice operators into the fractional derivatives 
with respect to coordinates. We represent these transformations by Figure 10. 
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From Lattice 


Qa,± 



to Continuum 

> 

ll } 



a Xj p 


TV 

/ \ 


Fourier series 


Inverse Fourier 

transform 


integral transform 



Figure 10: Transformation of lattice operators to fractional derivatives. 


Using the methods suggested in [60l[6T], we can prove the connection between the lattice 
operators and fractional derivatives of non-integer orders with respect to coordinates. 

The lattice operators fl2B]) . where K^{n — m) are dehned by flTTD and flTS]) . are transformed 
by the combination o him o into the fractional derivatives of order a with respect to 
coordinate Xi as 


^ o him o ( K 


a 

L i J 


d' 


iQ:,ib 


= a 


* d\xi 


(38) 


where = |aj| are the primitive lattice vectors, is the Riesz fractional derivative 

of order a > 0 with respect to Xi, and d°'~/d\xi\°' is the generalized conjugate Riesz derivative 
of order a > 0. The order of the partial derivative is defined by the order of lattice 

operator [“] and it can be integer and non-integer. 

Using the independence of the site vectors of lattice site ni = (ni,0,0), 112 = (0,^2, 0), 
ns = ( 0 , 0 , 77 . 3 ) and the statement fl38|) . we can prove that the continuum limits for the mixed 
lattice operators (HT} have the form 


o him o 


ai/3j 
i J 


* ^ d\xi\°‘^ d\xj\l^i' 


(39) 
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As a result, we obtain continuum limits for the lattice fractional derivatives in the form of the 
fractional derivatives of the Riesz type with respect to coordinates. 

The Riesz fractional derivative of the order a is dehned pm by the equation 


d\xi\^ di{m,a) 


1 


I I Q+l 
I \ 




(0 < a < m), 


(40) 


where (A™/)(zj) is a hnite difference of order m of a function /(r) with the vector step Zj = 
XiBi G for the point r G The non-centered difference is 


I f L I 

(Ar/)(z.)=X:(-l)‘^j£^/(r(41) 

and the centered difference 

m I 

(Ar/)(z,) = E(-l) h,(“lt)! -f*"- - - i) -.)■ (42) 


The constant di{m,a) is dehned by 

di{m, a) = 

where 




[a 


2"r(l + Q;/2)r((l + a)/2) sin(7rQ;/2) ’ 


7-1 


m\ 


AUa) = 2j2{-^y -u q, 

in the case of the non-centered difference flTTl) . and 




[m/ 2 ] I 

^m(a) = 2 . 1 , {'m/2 - 

in the case of the centered difference (IT2|) . The constants di{m,a) is different from zero for 
all Cl > 0 in the case of an even m and centered difference (A™/) (see Theorem 26.1 in 
P). In the case of a non-centered difference the constant di{m,a) vanishes if and only if 
a = 1, 3, 5,..., 2[m/2] — 1. Note that the integral (HD]) does not depend on the choice of m > a. 
The Fourier transform of the Riesz fractional derivative is given by 

Equation fH3|) can be considered as a dehnition of the Riesz fractional derivative of order a. 
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Using that (—= (—1)'^, the Riesz derivatives for even a = 2j are 




( 44 ) 


For a = 2 the Riesz derivative looks like the Laplace operator. The fractional derivatives 
/d\xi\'^ for even orders a are local operators. Note that the Riesz derivative /d\xi\^ 
cannot be considered as a derivative of first order with respect to |xj|. For a = 1 it is it looks 
like the square root of the Laplacian. The Riesz derivatives for odd orders a = 2j + 1 are 
non-local operators that cannot be considered as usual derivatives /dx^^^^. 

We also define the new fractional derivatives by the equation 


a"’- 

d\xi\°‘ 


dx, * 

A 

dxi 

d 


, dxi d\xi 


a—1 


0 < a < 1, 
a = 1, 
a > 1, 


(45) 


where d/dxi is the usual derivative of hrst order with respect to coordinate Xj, and “ is the 
Riesz potential of order (1 — a) (see Appendix) with respect to Xj, 


I-"“/(r) 


Ri-a{xi - Zi) /(r {Zi - Xi)ei)dzi, (0 < a < 1), 


(46) 


where e* is the basis of the Cartesian coordinate system. For 0 < a < 1 the operator (?“’“/c}|xj|“ 
is called the conjugate Riesz derivative [9]. Therefore, the operator 9"“/9|xi|" for all a > 0 
can be called the generalized conjugate Riesz derivative. 

The Fourier transform R of the fractional derivative fl45l) is given by 

^i ^aJP ) W = '*:.l*:ir''OT)(k) = »sgn(fc.)INV/)(k). (47) 

Using fl44D and fl4S]) . we get 


Q2j+i, ^ g^-^'+V(r) 


(48) 


The fractional derivatives (?“’“/c}|xi|" for odd orders a are local operators. Note that the 
generalized conjugate Riesz derivative d‘^~/d\xi\^ cannot be considered as a derivative of second 
order with respect to \xi\. The derivatives d^~/d\xi\°' for even orders a = 2j are non-local 
operators that cannot be considered as usual derivatives /dx^^. For a = 2 the generalized 
conjugate Riesz derivative is not the Laplacian. 
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Equations (Hill and fH8l) allow us to state that the usual local partial derivatives of integer 
orders are obtained from the operators in the following two cases: (1) for odd values 

a = 2j + 1 > 0 by d'^~/d\xi\'^ only; (2) for even values a = 2j > 0 by 8°"’^/d\xi\'^ only. The 
operators /d\xi\°' with integer odd a = 2j + 1 and d°'~/d\xi\°' with integer even a = 2j, 
where n G N, are non-local operators. Therefore we consider the lattice equations with the 


lattice operators IK and IK ’ 


as main lattice models to have the usual equations 


with local spatial derivatives in the case aj = /3j = 1 for all i = 1, 2, 3. 


4 Fractional Fokker-Planck equation for continuum 


Using the statements ([37]), fl38|) and fl3^ . where K~{n — m) are dehned by flTS]) . the lattice 
Fokker-Planck equation fl2^ are transformed by the combination o Limo into the 
fractional Fokker-Planck equation with derivatives of non-integer orders with respect to space 
coordinates. This space-fractional Fokker-Planck equation for the probability density f{r,t) 
has the form 


df{r,t) 

dt 


3 

E 

2=1 


D, 


a] 


Qoti, 

d\xi\^i 


/(r,t) 


3 3 


QoLi-i 

d\xi\°-i d\xj\^i 


(49) 


where Di{a) is the drift vector and Dij{a,P) is the diffusion tensor for the continuum that are 
dehned by the lattice coupling constants gi and gij by the relations 

Di{a) = af gi, Ai(a, /?) = 2 af gij. (50) 


Using the dehnition (051), the fractional Fokker-Planck equation (091) can be represented as 
the well-known continuity equation 


dt 


-E 


2=1 


dJi{r,t) 

dxi 


where J* is the probability how 


i=i 


Ji{r,t) = < 


d\xj\f^^ 
8^^ 


Di{a)f{r,t) - -J2Dij{a,(^) ift. 

j=i J 


8\xj 




8\xj\^^~^ 8\xj\l^^' 


( 51 ) 


1 ^ 

A(a)i;-"7(r.*) - 5 $^ A7a./l)iy"‘75nr^/(r7) 0 < a. < 1, 


CXi 1, 


piOii 1 )+ f (Y 1 - ^ ^ 

- 7;^Dij{a,P) «* > 1 , 


( 52 ) 


14 















Note that coincidence of orders of fractional derivatives in the hrst and second terms allows 
us to represent the fractional Fokker-Planck equation in the form of the space-fractional 
continuity equation. The fractional Fokker-Planck equation fH9|) can be represented as the 
fractional continuity equation 


dt 

where is the probability flow 

Afrac) 


E 

i=l 


Qai, 

d\xi\°‘i 


(53) 


J 


dl^i 


(r,() = Di(a)f(r,t) - 


(54) 


If Qfj = 1, the continuity equation fl5^ has the standard form. 

For one-dimensional case with Di(a) = 0 and /(r,f) = f{x,t), equation can be repre¬ 
sented in the form 

= K(,i)Wf(x,t), (55) 

where K{fi) is the generalized diffusion constant, 


*’( a ) = 2^11(0,^), 


(56) 


and is the fractional derivative of order /i, 




^ai,- g0i,- 

"oi I "HI I fl, ) h CXi -f- Pi. 

a|a:|"i o \ x \ p ^ 


(57) 


Note that for sufficiently good functions, the operator (IFr|) can be represented in the form 
V/4 — d^'~^/d\x\^, but it cannot be done in the general case. Equation fl5^ describes the 
fractional diffusion processes with the Poissonian waiting time and the Levy distribution for 
the jump length (see Section 3.5 of [1]). In [1] the space-fractional diffusion equation fl55ll 
contains the Weyl fractional derivative of order /r, which is equivalent to the Riesz operator 
d^’^/d\x\^ in one dimension. The solution of equation fl55|l can be obtained analytically by 
using the Fox function H 22 (for details see Section 3.5 in |1] and PI). The exact calculation 
of fractional moments [1] gives 




2(iF(/i))^//-F(-V/i)F(l + ^) 

pF(-5/2)F(l + 5/2) 


where 0 < 5 < p < 2. 


( 58 ) 
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The time-fractional lattice Fokker-Planck equation fl30l) are transformed by the combination 
o Limo into the space-time fractional Fokker-Planck equation 


dt 


A(a) 


2=1 


QOLi 1 

d\xi\°‘^ 






Qai , Ql3j , 


RL ni-7 
0 




(59) 

where is the Riemann-Liouville fractional derivative with respect to time that describes 

the power-law memory. For one-dimensional case with Di{a) = 0 and /(r, t) = f{x, t), equation 
(15^ can be represented in the form 


(60) 

where 

^(h,7) = (61) 

and the fractional derivative is dehned by fl57)) . Equation fl60|) describes a random walk 
characterized by waiting time and jump length (see Section 3.6 in [1]). The competition between 
long rests (waiting events) and long jumps (motion events) in the Levy walks processes is given 
[71] as 

( 0 < 7 < 1 , 

{x^{t)) ~ (62) 

[ 7 > 1, 

where 1 < /i < 2. It should be noted that the continuum form of the Levy flights is described 
by the drift term with the hrst order derivative (oj = 1) as proposed in [39] and derived from 
the continuous time random walk in [36]. The solutions of Cauchy problems for the space-time 
fractional diffusion equation with the Riesz-Feller fractional derivatives are described in [Tl] . 
The time-fractional lattice Fokker-Planck equation fl3T|) is transformed by the combination 
o Lim o into the space-time fractional continuum Fokker-Planck equation 

(63) 

where 

^CFP = ° Lim o (64) 

is the continuum Fokker-Planck operatorof the form 


ra,/? _ 

^CFP ~ 


J2d,(c) 


1 


2 = 1 


d\xi 


1 3 3 OQ,, 




i=l j=l 


d\xi\^^ d\xj\l^i 


(65) 


For CKj = /3j = 1, equation fl63|) takes the form of the time-fractional Fokker-Planck equation 
that is suggested in [361 ETJ [38| . 
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If a* = = 7 = 1 for all z = 1 , 2 , 3, then equations (09]), fl5^ and fl 6 ^ for the probability 

density f{r,t) give the well-known Fokker-Planck equation in the form 


dt 


3 „ 3 3 




dxi dx. 




( 66 ) 


i=l “ ‘ i=l j=l ‘ ^ 

where Di = Di{l) is the usual drift vector and Dij = -Dp(l, 1 ) is the usual diffusion tensor. 


5 Conclusion 

In this paper three-dimensional lattice models with long-range drift and diffusion of particles 
are suggested. These proposed lattice models can be considered as a possible microscopic basis 
to describe the anomalous diffusion in continuum. The suggested type of lattice long-range 
drift and diffusion can be considered for non-integer (fractional) values of the parameters 
/3j, 7 . This allows us to have lattice equations for the fractional nonlocal diffusion and transport 
processes. The proposed forms of the drift and diffusion of particles in lattice allow us to obtain 
the continuum equations with the generalized conjugate Riesz derivatives of fractional orders 
by using the approaches and methods proposed in [ 6 OI El]. The suggested three-dimensional 
models with long-range lattice drift and diffusion of the types flT7|) and fflSj) can be considered as 
lattice analogs of the fractional diffusion and drift in nonlocal continuum. Different fractional 
generalizations of Fokker-Planck equation for continuum can be obtained by using the sug¬ 
gested lattice approach. We expect that the proposed three-dimensional lattice Fokker-Planck 
equations can play an important role in the description of nonlocal processes in microscale and 
nanoscale because at these scales the interatomic interactions can be prevalent in determining 
the properties of media. 

Let us note some possible generalizations of the proposed lattice models. We assume that the 
suggested lattice Fokker-Planck equations can be generalized in the form of lattice Kramers- 
Moyal equation for the case of the high-order terms by using the different fractional-order 
derivatives. The suggested lattice models can be generalized for the lattices with dislocation 
and disclinations that are connected with non-commutativity of the lattice operators fl26|l . In 
this paper, we consider the primitive orthorhombic Bravais lattice for simplihcation. It is 
interesting to generalize the suggested consideration for other type of Bravais lattices such as 
triclinic, monoclinic, rhombohedral and hexagonal. The suggested models of unbounded lattices 
can be generalized for the bounded physical lattices. We also assume that the proposed lattice 
approach to the fractional diffusion can be generalized for lattices, which are characterized by 
fractal spatial dispersion 1751 [761 [77], and correspondent models for fractal media |78l [79] (see 
also [801EI1182I). 

We also can note some remaining challenges and questions in the suggested approach to 
fractional diffusion. The function f{n,t) has the meaning of probability density on the lattice, 
and it should be positively dehned. It is known that this condition for continuum case leads 
to restriction for the parameters a, f3, 7 , fi. For example, we have the condition 0 < /i < 2 
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for Levy processes on continuum. A rigorous consideration of positiveness of /(n, t) for set of 
these parameter does not exist at the present time. Exact mathematical conditions of existence 
of solutions for the lattice Fokker-Planck equation can be important to describe anomalous 
long-range particle drift and diffusion on three-dimensional physical lattices. 
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Appendix: Riesz fractional integral 

The Riesz fractional integration is dehned by 

(67) 

The fractional integration (l67|) can be realized in the form of the Riesz potential dehned as the 
Fourier’s convolution of the form 

^xfi^)=[ Ra{x - z)f{z)dz, (a > 0), (68) 


RJx) = 


(69) 


where the function Ra{x) is the Riesz kernel. If a > 0, the function Ra{x) is dehned by 

7 “^(a)|a;|"“” a^n + 2k, 

— 7 “^(q:)|x|"“” In |x| a = n + 2k, 
where n G N, and the constant 7n(Q:) has the form 

2“7r”/2p^^/2)/P(^) a^n + 2k, 

(_l)D-«)/ 22 «-i 7 rV 2 r(a/2) P(1 + [a - n]/2) a = n + 2k. 

The Fourier transform of the Riesz fractional integration is given by 

R(rj{x)) = \knRf){k). (71) 


7n(a) = 


(70) 
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Figure 1: Plot of the function F+(a:, y) ffT^ for the range x G [0, 6] and y = a E [0, 6]. 



Figure 2: Plot of the function F_{x, y) (j2nil for the range x G [0, 6] and y = a E [0, 6]. 
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Figure 3: Plot of the function F+(x, y) flT^ for the range x G [5,15] and y = a E [0, 6]. 



Figure 4: Plot of the function F_{x, y) (j2nil for the range x G [5,15] and y = a E [0, 6]. 
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Figure 5: Plot of the function F+(a:, y) ffT^ for the range x G [1, 5] and y = a ^ [1)4], 



Figure 6: Plot of the function F_{x, y) (j2nil for the range x G [1, 5] and y = a E [1,4]. 
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Figure 7: Plot of the function F^{x) (IT^ with a = 1.5 and a = 1 for the range x G [0, 7]. 



Figure 8: Plot of the function F-{x) (120]) with a = 1.5 and a = 1 for the range x G [0, 7]. 
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